arXiv: 1504.08243vl [math.DS] 30 Apr 2015 


Exponentially small asymptotic formulas for the length spectrum in some billiard tables^ 


P. Martin'’, R. Ramfrez-Ros^ * *, A. Tamarit-SarioP 


^Departament de Matemdtica ApUcada I, Universitat Politecnica de Catalunya, Diagonal 647, 08028 Barcelona, Spain 
^Departament de Matemdtica ApUcada IV, Universitat Politecnica de Catalunya, Campus Diagonal Nord, Edifici C3. C. Jordi Girona, 1-3, 08034 Barcelona, Spain 


Abstract 

Let ^ > 3 be a period. There are at least two (1, ^)-periodic trajectories inside any smooth strictly convex billiard table, and all 
of them have the same length when the table is an ellipse or a circle. We quantify the chaotic dynamics of axisymmetric billiard 
tables close to their borders by studying the asymptotic behavior of the differences of the lengths of their axisymmetric (1,^)- 
periodic trajectories as o' —> +oo. Based on numerical experiments, we conjecture that, if the billiard table is a generic axisymmetric 
analytic strictly convex curve, then these differences behave asymptotically like an exponentially small factor q'^^e times either a 
constant or an oscillating function, and the exponent r is half of the radius of convergence of the Borel transform of the well-known 
asymptotic series for the lengths of the (1, q')-periodic trajectories. Our experiments are restricted to some perturbed ellipses and 
circles, which allows us to compare the numerical results with some analytical predictions obtained by Melnikov methods and 
also to detect some non-generic behaviors due to the presence of extra symmetries. Our computations require a multiple-precision 
arithmetic and have been programmed in PARJ/GP. 
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1. Introduction 

Billiards as a dynamical system go back to Birkhoff lU]. 
Let g be a closed smooth strictly convex curve in the Euclidean 
plane. The Birkhoff billiard models the motion of a particle in¬ 
side the region enclosed by Q. The particle moves with unit 
velocity and without friction following a straight line; it reflects 
elastically when it hits Q. Therefore, billiard trajectories con¬ 
sist of polygonal lines inscribed in Q whose consecutive sides 
obey to the rule “the angle of reflection is equal to the angle 
of incidence.” Such trajectories are sometimes called broken 
geodesics. See for a general description. 

A {p, ^)-periodic billiard trajectory forms a close^olygon 
of q sides that makes p turns before closing. Birkhoff yj] proved 
that there are at least two different Birkhoff {p, q')-periodic bil¬ 
liard trajectories inside Q for any relatively prime integers p 
and q such that 1 < p < ^. 

The length spectrum of Q is the subset of K.+ defined as 
£S{Q) = /N U IJ 

(P^4) 

where / = Length(2) and c IR+ is the set of the lengths 
of all (p, ^)-periodic billiard trajectories inside Q. The maximal 
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difference among lengths of (p, ^)-periodic trajectories is the 
non-negative quantity 

ACa?) =supA<^’«)-infA'^’'?'. 


Many geometric and dynamical properties are encoded in 
the length spectrum nS{Q) and the differences 

An old geometric question is: Does the set -LS{Q) allow 
one to reconstruct the convex curve Q? The length spectrum 
and the Laplacian spectrum with Dirichlet boundary conditions 
are closely related 13]. Therefore, the question above can be 
colorfully restated asldj: Can one hear the shape of a drum? 
We refer to the book ITt] for some results on this question. 

The difference is important from a dynamical point 
of view, because it is an upper bound of Mather’s AWp/g. In 
its turn, AWp/g is equal to the flux through the (p, q')-resonance 
of the corresponding billiard map fS-llj. Thus, the variation 
of A^P’‘i'> in terms of the rotation number pjq 6 (0,1) gives in¬ 
formation about the size of the different chaotic zones of the 
billiard map. See Section |2] for a more complete description of 
these ideas. 

Here, our main goal is to gain some insight into the billiard 
dynamics close to the boundary of the billiard table. We focus 
on the (1, ^)-periodic billiards trajectories; that is, we set p = 1. 
We want to determine the asymptotic behavior of 


A"’«) = supV’«^-infA"’«) 


q ^ H-oo. 

Let e be the length of a (1, ^)-periodic billiard 
trajectory inside Q. It does not matter which one. Marvizi and 
Melrose lll2ll proved that if Q is smooth and strictly convex. 
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then there exist some asymptotic coefficients /, = Ij(Q) such 
that 

X ^ oa. (1) 

J>0 


for some jj. e {0,1), some amplitude a + 0, and some phases 
(fj e R, when these singularities are 

+di (if and only if /i = 1), +ySi ± 6i,... ,±l3j ± 5i. 


For instance, Iq - I - Length(2) and /i = , 

where k and ds are the curvature and the length element of Q, 
respectively. The symbol x means that the series in the right 
hand side is asymptotic to The asymptotic coefficients Ij 
can be explicitly written as integrals over Q of suitable algebraic 
expressions of k and its derivatives. The first five coefficients 
can be found in Q. The asymptotic series O does not depend 
on the choice 6 so 

lim =0, 'ik> 0. 

»+oo 

That is, the differences are beyond all orders in q. In fact, 
they satisfy the following exponentially small upper bound in 
the analytic case ill. If 2 is a closed analytic strictly convex 
curve, then there exist constants K, qo,a > 0 such that 

^(1.9) < Mq > qo. 


The exponent a is related to the width of the complex strip 
where a certain 1-periodic angular coordinate is analytic. If 
a billiard map (or any analytic exact twist map) has a rotational 
invariant circle of Diophantine rotation number w, there exist 
other exponentially small upper bounds for (or for the 
residues of (p, ^)-periodic orbits) when pjq —» u. See 01514171. 

Similar singular behaviors have been observed in problems 


about the splitting of separatrices of analytic maps 1118142911 . All 


these splittings are not exponentially small in a discrete big pa¬ 
rameter ^ e N, but in a continuous small parameter h > 0. 
Namely, h is the characteristic exponent of the hyperbolic fixed 
point whose separatrices split. Thus, we may think that h - Ijq 
for comparison purposes. The splitting size in many analytic 
maps satisfies the exponentially small asymptotic formula 


'‘splitting size” “ A(l//r)/r “e 


0 ^ ( 2 ) 


for some exponent r > 0, some power m 6 R, and some func¬ 
tion A( 1 /h) that is either constant or oscillating. The exponent r 
and the function A(l// 2 ) are determined by looking at the com¬ 
plex singularities closest to the real axis of the homoclinic so¬ 
lution of a limit Hamiltonian flow related to the map. Such 
methodology has been rigorously established for the standard 
map 0221] ■ the Henon map 12411 . and some perturbed McMil¬ 
lan maps 1^ ^ 3- It has also been numerically checked 
in certain billiard maps and several polynomial maps 1^ . 
but there are other maps where it fails 02911 . Let us briefly re¬ 
call some claims about polynomial standard maps contained 
in MM- First, r - 2n5, where d is the distance of these 
singularities to the real axis. Besides, 


j 

A(l//!) = pa 12 + a ^ cos(27r/3y//! + ipj), (3) 

1=1 


For instance, the limit Hamiltonian flow for the standard map 
is a pendulum, so +7ri/2 are the closest singularities to the real 
axis and the “splitting size” is the so-called Lazutkin constant 
0)0 ^ 1118.827706 times see B . 

It is also known that, usually, r - p/2, where p is the radius 
of convergence of the Borel transform of the di vergent asymp¬ 
totic series that approaches the separatrices 1 2j. , 24 - 2^ . 

By looking at our billiard problem from the perspective of 
those results (and others not mentioned here for the sake of 
brevity), it is natural to make the following conjecture. This 
conjecture is strongly supported by our numerical experiments. 


Conjecture 1 . If Q is a closed analytic strictly convex curve, 
but it is neither a circle nor an ellipse, the asymptotic series (O} 
diverges for all period q' 6 N, but it is Gevrey-1, so its Borel 
transform 






AO 


( 2 /- 1 )!’ 


(4) 


has a radius of convergence p 6 (0, -too). Set r — p/2. 

If Q is a generic axisymmetric algebraic curve, then 


X \A{q)\q-^t-'‘‘, q ^ -Hoo, (5) 


for some function A{q) that is either constant: A{q) — a/2 ^ 0, 
or oscillatory: 


j 

A{q) - pal2 + a'^ cos(27Tfjq), (6) 

i=i 

with p 6 (0,1), a 0, 7 > 1, and 0 < < ■ ■ ■ < fj. The cases 

A(q) — a/2 andA{q) — a cos(2nf3) take place in open sets of the 
space of axisymmetric algebraic curves. All the other cases are 
phenomena of co-dimension one. 

If Q is a generic bi-axisymmetric algebraic curve, A^’t^ has 
the previous asymptotic behavior when q is even and q —> H-oo, 
but A*'’?) = 0{q-\-^''‘>) when q is odd and q —> -Hoo. 

We stress that the oscillating function Q has some phases, 
but there are no phases in the oscillating function This 
phenomenon is not new. The asymptotic formulas for the ex¬ 
ponentially small splittings of generalized standard maps with 
trigonometric polynomials do not have phases either 12611 . 

A curve is axisymmetric when it is symmetric with respect 
to a line, and bi-axisymmetric when it is symmetric with respect 
to two perpendicular lines. A planar curve is algebraic when its 
points are the zeros of some polynomial in two variables. We 
require strict convexity, since it is already an essential hypoth¬ 
esis in the smooth setup. We only consider algebraic curves by 
comparison with the above results about polynomial standard 
maps. Our algebraic curves have no singular points, because 
we ask them to be closed and analytic. 
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If 2 is a circle of radius ro, all its (p, ^)-periodic billiard 
trajectories have length 2roqsm{Trpfq), so = 0 for all 

piq e (0,1), and the asymptotic series ([TJ becomes 

= 2roq sinin/q) = 2ro V q~^^, 


which converges for all q. In particular, p - +oo. Ellipses have 
analogous properties. This has to do with the fact that elliptic 
and circular billiards are integrable. A conjecture attributed to 
Birkhoff claims that the only integrable smooth convex billiard 
tables are ellipses and circles 1^ . Following the discussion on 
the Mather ^S-function contained in lll3n . this old conjecture is 
reformulated as; The series in m converges for some period 
O' 6 N if and only if Q is an ellipse or a circle. 

In this paper, we present several numerical experiments and 
some analytical results that support Conjecture [T] For the sake 
of simplicity, all numerical experiments are carried out using 
the model tables 


information when n is odd and q even; = 0 in such case. 
See Propositionl^for details. 

Which is the relation between the asymptotic formula Q 
and the first order Melnikov computation (|8]l? The answer is 
that r + c and m„ + —3, so the Melnikov method does not 
accurately predict the singular behavior Nevertheless, 

lim^^o r = c, so some information can be retrieved from the 
Melnikov method, at least for perturbed ellipses. 

The case of perturbed circles is harder. See Section 

Symmetries are another reason for the choice of tables (I7]i. 
On the one hand, symmetries greatly simplify the computa¬ 
tion of periodic trajectories. To be precise, we just compute 
the signed difference Dq between two particular axisymmetric 
(1, q')-periodic trajectories, instead of or AWi/^. Clearly, 
\Dq\ < Often, \Dq\ - = AWi/^. See Proposition |2] 

On the other hand, bi-axisymmetric curves are a very particular 
class of axisymmetric curves, so our model tables may display 
other asymptotic behaviors when n is even. We will check that 
this expectation is fulfilled. Concretely, 


Q^{{x,y)eM.^:x^+y^/b^ + ey’'^l}. (7) 


A''’-?) X \B(q)\q-^e-^^f q ^ -Hoo, 


Here, e (0,1] is the semi-minor axis, e 6 R is the perturbative 
parameter, and n 6 N, with 3 < « < 8, is the degree of the 
perturbation. We will refer to Q as a perturbed ellipse when 
Q < b < \ and as a perturbed circle when b - \. Next, we 
explain the four main reasons for this choice of billiard tables. 

As a first reason, we know that all the billiard tables dill are 
nonintegrable for n > 3 and 0 < e <sc 1, and so the dynamics 
inside them should be far from trivial. The question of which 
perturbed ellipses give rise to integrable billiards is addressed 
in 03111 . Theorem 4.1 of that paper imply that the tables (|7]i are 
nonintegrable if0<f7< l,n>4is even, and e is small enough. 
This result can be extended, after some technicalities, to odd 
degrees. Furthermore, all integrable deformations of ellipses 
of small eccentricities —this includes, of course, circles— are 
ellipses so the tables (|7]i are nonintegrable if b - 1, « > 3, 
and e is small enough. 

The second reason is that we want to use some Melnikov 
methods that are well suited for the study of billiards inside 
perturbed ellipses and perturbed circles 03311341] . We recall that 
A(1.?) - 0 for any q' > 3 in elliptic billiards. Thus, since the 
difference = A^'’^'(e) is analytic in e and vanishes at e = 0, 
we know that 


for some coefficient Aj*’'*' e R that can be computed explicitly. 
To be precise, it turns out that if 0 < /? < 1 then 

Aj*’'** M„q'"''e^‘^'^, q -too, (8) 

for some Melnikov exponent c > 0 not depending on n, some 
Melnikov power m„ 6 Z, and some Melnikov constant M„ ^ 0. 
These three Melnikov quantities can be explicitly computed, 
but we have carried out the computations only for the cubic 
(n = 3) and quartic (n = 4) perturbations for the sake of brevity. 
Besides, lim^^i c - -too. The Melnikov method provides no 


for some constant or oscillating function B(q) when n is even 
and q is odd. This asymptotic behavior has several differences 
with respect to the generic one conjectured in Q. Both the 
exponent in e and (if any) the frequencies 0 < fi < ■ ■ ■ < /3j 
are doubled, the power in q~^ is increased by one, etcetera. We 
think that this new asymptotic behavior is generic among bi- 
axisymmetric algebraic curves when the period q is odd. 

The last reason for the choice of such simple billiard tables 
is to reduce the computational effort as much as possible. In 
particular, we limit the degree of the perturbation to the range 
3 < « < 8 for this reason. Recall that each set is con¬ 
tained in an exponentially small (in q) interval, so the compu¬ 
tation of (or Dq) gives rise to very strong cancellations. 
This forces us to use a multiple-precision arithmetic to com¬ 
pute them. We have performed some computations with more 
than twelve thousand digits, based on the open source PARI/GP 
system ll35n . Similar computations in the setting of splitting of 


separatrices of analytic maps can be found in 112IL |25 


Finally, we recall that the exponent r is found by looking at 
the complex singularities of the homoclinic solution of a limit 
Hamiltonian flow in many cases of splitting of separatrices. 
Does such kind of limit problem exist in our billiard setting? 
Unfortunately, we do not have a completely satisfactory answer 
yet, but we propose a candidate in Section |3] It is empirically 
derived by using the Taylor expansions of the billiard dynam¬ 
ics close to the border given by Fazutkin in ll^ . Fet k{s) be 
the curvature of Q in some arc-length parameter s 6 R//Z. Fet 
f e R/Z be a new angular variable defined by 




C = 




,2/3 


di. 


(9) 


Fet 6 be the distance of the set of singularities and zeros of the 
curvature K(f') to the real axis. We thought that r - 2n6, but our 
experiments disprove it. We have only obtained that r < 2n6, 
the equality being an infrequent situation. But there are some 
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good news about our candidate. First, the Melnikov exponent 
is c = 27t6, when Q - {(.r,y) e +y^lb^ = 1), with 

0 < b < 1. See Proposition H] Second, we have also seen that 
if b = 1 and « > 3 is fixed, then there exist some constants 
6 R,< 77„, such that 


r — 


I log el 

n 


+;r« + 0(1), 


2n5 - 


I log el 

n 


+ rjn + o(l). 


as e ^ 0^. The second formula is proved in Proposition |7] the 
first one is numerically checked in Section |6l Therefore, our 
candidate captures exactly the logarithmic growth of the expo¬ 
nent r for perturbed circles. Third, our experiments suggest that 
r = 2n6 when b - 1, n e {5,7), and e 6 (0,1/10). 

The paper has the following structure. Section |2 contains 
the dynamical interpretation of Mather’s AWp/q. We discuss 
our candidate for limit problem in Section[3 The axisymmetric 
tables and their axisymmetric periodic billiard trajectories are 
presented in more detail in Section |4] The main results about 
perturbed ellipses and perturbed circles are described in Sec¬ 
tions |5] and |6] respectively. All proofs have been relegated to 
the appendices. 


2. Twist maps, actions, Mather’s AW, and billiards 

We recall some results about exact twist maps and billiards. 
We refer to the books iM] and the surveys iflH[Til] for a more 
detailed exposition. 

Let T = R/Z and I = (j_,y+) c R for some -oo < y_ < 
y+ < +00. Let oj - dx A dy be the canonical area form on T x /. 
Note that u - — dA, where A - y dx. A smooth diffeomorphism 
/:Tx/^Tx/isan exact twist map when it preserves a>, 
has zero flux, and satisfies the twist condition ^ > 0, where 
F:Rx/—»Rx/, F(x,y) — (xi,yi), is a (fixed) lift of /. 

We also assume that / can be extended as rigid rotations 
of angles g± to the boundaries C± - T x {y±). We know that 
Q- < g+ from the twist condition. Let p and q be two relatively 
prime integers such that pjq e (p-,p+) and q > \. A point 
(x,y) 6 R X / is {p, q)-periodic when F‘^{x,y) = (x -H p,y). The 
corresponding point (x, y) e T x / is a ^'-periodic point of / that 
is translated p units in the base by the lift. A {p, q')-periodic 
orbit is Birkhojf when it is ordered around the cylinder in the 
same way that the orbits of the rigid rotation of angle pjq. The 
Poincare-Birkhoff Theorem states that there exist at least two 
different Birkhoff {p, ^)-periodic orbits llliS]. 

Let E - {(x, xi) € R^ : < xi - x < p+). Then there exists 

a function : £ —» R such that hix -H 1, xi -H 1) = h{x, xi) and 


yi dxi — y dx = d/ 2 (x, xi). 



Note that, although the series for W[x] may be divergent, 
only involves two terms of the series, and so VW is well defined. 
\f O - {{xj,yj)]jez is a {p, ^)-periodic orbit of /, then 

h{Xj+q,Xj+q+\) = h{Xj + p,Xj+\ + p) ^ h{Xj,Xj+\), 

so there are only q different terms in the action functional W, 
which encode the (p, q')-periodic dynamics. In particular, any 
ip, ^)-periodic orbit O = {{xj,yj)}jei, is in correspondence with 
a stationary configuration x = (xq, ... , x^_i) e R^“^ of the 
(p, q)-periodic action 

^(P,9)[x] = /2(xo,Xi) -t- /2(Xi,X2) H - + h{Xq^i,Xo + p). (10) 

We say that is the (p, ^)-periodic action of 

the (p, g)-periodic orbit O. The Birkhoff (p, ^)-periodic orbits 
provided by the Poincare-Birkhoff Theorem correspond to the 
minimizing and minimax stationary configurations of . 

Mather defined the quantity AWpjq > 0 as the action of 
the minimax periodic orbit minus the action of the minimiz¬ 
ing one 1^. Mather’s AWpjq has a nice dynamical interpre¬ 
tation. It is equal to the flux through any homotopically non 
trivial curve without self-intersections passing through all the 
points of both the minimizing and the minimax (p, q')-periodic 
orbits. This is fhe. MacKay-Meiss-Percival action principle |@]. 
Thus, AWpjq gives a rough estimation of the size of the (p, q)- 
resonance of the twist map. We also recall that the hyperbolic 
(respectively, elliptic) periodic orbits in a given resonance are 
generically minimizing (respectively, minimax). 

Next, we adapt these ideas to billiard maps. 

Let 2 be a smooth strictly convex curve in the Euclidean 
plane. For simplicity, we assume that I = Length(2) = 1. 
Let 7 : T —> 2, i 7('S), be an arc-length counterclockwise 
parametrization of Q. The bounce position of the particle inside 
Q is determined by the arc-length parameter s. The direction of 
motion is measured by the angle of incidence r e (0, n). Let 


This function is called Lagrangian or generating function. It 
is determined modulo an additive constant. Twist maps satisfy 
the following classical Lagrangian formulation. Their orbits 
are in one-to-one correspondence with the (formal) stationary 
configurations of the action functional 

R* 3 X = {Xj)jeZ I-+ W[x] = ^ h{Xj,Xj+i). 

jeZ 


/:Tx(0,7r)^Tx(0,7r), f(s,r)^(suri), (11) 

be the corresponding billiard map. Figure[T]illustrates this map. 
The coordinates (s, r) are called Birkhoff coordinates. 

Let i{s, si) = l 7 (s) - 7 (si)| be the Euclidean distance be¬ 
tween two impact points on Q. It is easy to prove that 


4 


— (i, ii) = -cosr, 
os 


—(i, ii) = cosri. 
ds\ 








ip,q) 



^^PKi 

( 1 , 2 ) 

4.000000 

3.828482 

0.171577 

(1,4) and (3,4) 

5.594652 

5.536901 

0.057751 

(1,3) and (2,3) 

5.115169 

5.112940 

0.002229 

(3,8) and (5,8) 

14.773311 

14.772302 

0.001009 

(5,12) and (7,12) 

23.151909 

23.150969 

0.000940 

(3,10) and (7,10) 

15.925337 

15.924445 

0.000892 

(1,6) and (5,6) 

5.904338 

5.903527 

0.000811 

(2,5) and (3,5) 

9.366997 

9.366503 

0.000494 

(1,8) and (7,8) 

6.024507 

6.024232 

0.000275 

(3,7) and (4,7) 

13.455442 

13.455236 

0.000206 

(1,5) and (4,5) 

5.785133 

5.785011 

0.000122 


Table 1: The biggest Mather’s for the billiard inside + 7“*/10 = 1. 

If we consider the coordinates (x, y) - {s, - cos r) e T x (-1,1), 
then the billiard map / becomes a smooth exact twist map with 
Lagrangian !i{x, xi) = -^(s, si) and boundary rotation numbers 
= 0 and = 1. That is, the action of a periodic billiard 
trajectory is, up to the sign, its length. In particular, Mather’s 
is the length of the (p, q')-periodic billiard trajectory that 
minimizes the action (and so, maximizes the length) minus the 
length of the minimax one. Generically, 

AWp/g = 

where and are the lengths of the hyperbolic and el¬ 
liptic (p, ^)-periodic billiard trajectories inside Q. For instance, 
= 4a, Zid.Z) _ 4 ^^ aWi /2 - 4(a - b) for the billiard 
inside the ellipse x^/a^ - 1 - jb^ - 1 with 0 < b < a. 

Note that any (p, ^)-periodic billiard trajectory gives rise to 
a iq - p, ^)-periodic one by inverting the direction of motion. 
This means that AWp/q - AW(g-p)/g for all pjq e (0,1 /2). 

We have listed the biggest Mather’s AWp/^ for the billiard 
inside the perturbed circle + y^ + y‘^110 - 1. See Table [T] 
The rest of Mather’s AWp/^ are smaller that lO^”*. The values 
in the table suggest that the (1,2)-resonance and both (p,4)- 
resonances should be the most important ones. This prediction 
is confirmed in Figure |2] where we display the biggest reso¬ 
nances of the billiard map inside x^ H- -H y ^/10 = 1 . 

Mather’s ISM pjq allow us to single out the most important 
resonances, but they do not give an exact measure of the size of 
resonances. To begin with, there is not a unique way to define 
such size. A choice is the area Apjq of the Birkhoff instability 
region that contains the (p, q')-resonance. A Birkhoff instability 
region is a region of the phase space delimited by two rotational 
invariant curves (RICs) without any other RIC in its interior. If 
we have a twist map with a (p, q')-resonant RIC, then AW pjq - 
0(e) and A„in - under generic perturbations of order 

0(e). See 11371] . This shows up a clear difference between these 
two quantities. For instance, the billiard map inside the circle 
x" -I- y^ = 1 has a (1,2)-resonant RIC, which is destroyed under 
the perturbation x^ -H y^/(l - = 1. However, this perturbed 

billiard table is integrable (it is an ellipse), so both quantities 
can be analytically computed; AIT 1/2 = 4e and A 1/2 = 

We omit the details. 



Figure 2: The biggest (p,(j)-resonances of the billiard map f{s,r) = (si,ri) 
inside the perturbed circle Q = {(jr, )^) € + v'^/lO = 1}. We recall 

that / = Length(Q). All (p, ^)-resonances with odd period q have 2q elliptic 
islands due to the bi-axisymmetric character of the curve. From bottom to top: 
(1,8), (1,6), (1,5), (1,4), (3,10), (1,3), (3,8), (2,5), (5,12), (3,7), (1,2), and 
their {q - p, q) symmetric counterparts. 

3. A candidate for limit problem 

To begin with, we recall how to obtain the limit problem 
for the splitting of separatrices of the generalized standard map 
/(x,y) = (xi,yi) given by 

xi=x-l-yi, yi = y -I- ep(x). (12) 

For simplicity, we assume that p{x) is a polynomial, p(0) = 0, 
and p'( 0 ) = 1 , so the origin is a hyperbolic fixed point of / 
with eigenvalues A = e^’ and 4“' = where e = 4 sinh^(/r/2). 
There is numerical evidence that the splitting size in this kind of 
polynomial standard maps satisfies the asymptotic formula (| 2 ]i 
for some exponent r > 0, some power m e R, and some con¬ 
stant or oscillating function A( 1 / h). We determine the exponent 
following lEill . 

First, we transform the original map into the map 

xi-x + pzi, z\-z+pp{z) 

by means of the scaling z- yip, where p - ffe. Note that p ^ h 
as e ^ 0^. The dynamics of this map for small p resembles 
the dynamics of the p-time flow of the Hamiltonian Ho(x, z) - 
z^/2-J p(x) dx. Besides, the origin is a hyperbolic equilibrium 
point of the Hamiltonian system 

x' = d^Ho(x,z) = z, y' = -dMo(x,z) = p{x). 

If the singular level set {(x, z) 6 : Hq(x,z) = Ho(0,0)} con¬ 

tains a separatrix to the origin, then we compute the flow on it 
and we get a homoclinic solution (xo(f), Zo(^)) that can be seen 
as the limit of the map on its separatrices when e ^ 0^. Such 
homoclinic solution is determined, up to a constant time shift, 
by imposing 

= PMO), lim xo(^) = 0 . 

^—>±co 
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It turns out that there exists 5 > 0 such that xo(f) is analytic 
in the open complex strip = {^ 6 C : |3^| < 5) and has 
singularities on the boundary of Is- Then r - 2ji6. This claim 
is contained in ll23l] . although a complete proof is still pending. 
However, Fontich and Simo proved the following exponentially 
small upper bound in If Q- e (0, S), then there exist some 
constants K, ho > 0 such that 

“splitting size” < VO < /j < ho- 

We want to emphasize an essential, but sometimes forgotten, 
hypothesis of the Fontich-Simo Theorem. Let 

o-oi^) = (vo(f),yo(f)) 

be the original homoclinic solution. The generalized standard 
map should have an analytic extension to a complex neigh¬ 
borhood in of (Toila)- If p{x) is a polynomial, then / can 
be extended to the whole and this hypothesis is automati¬ 
cally fulfilled. On the contrary, it remains to be checked when 
p : K —> K is just a real analytic function. 

Next, we adapt these ideas to our billiard problem. 

Let Q be an analytic strictly convex curve in the Euclidean 
plane. Set I - Length(2). Let k{s) be the curvature of Q in 
some arc-length parameter s e R/IZ. Note that k{s) > 0 for 
all s e R/IZ. Let p(s) - \Ik{s) be the radius of curvature. We 
are interested in the dynamics of the billiard map (fTTT) when the 
angle of incidence r tends to zero. More precisely, we consider 
that r - 0(\lq) and q —> -i-oo. 

Lazutkin 13611 gave the Taylor expansion 

Si = s -I- 2Q{s)r + 0(r“), n - r- 2p'(s)r^/3 + O(r^) 

for the dynamics of the billiard map (fTTT i around r = 0. Once 
fixed a period q » 1, we take p - 1/q' <s 1 as the small 
parameter. Then we transform the previous expansion into 

2 

Si = s -I- pq{s)v^'^ + Oit?), vi = V - -pq{s)v^^^ + Oij?), 

by means of the change of variables ^Jv - 2rfp. The billiard 
dynamics for small p resembles the dynamics of the /r-flow of 
the Hamiltonian Ho(s, v) = That is, the yu-flow of the 

Hamiltonian system 

s' = e(s)v'^^ v' = 

We compute the flow on the level set Kq {Ho{s, v') s |C^), 
for some constant C > 0. If (s, v) e Kq, then the first equation 
of the Hamiltonian system reads as 


approximated the billiard dynamics by the yu-time of the Hamil¬ 
tonian flow with p - Ijq. Any (1, q')-periodic trajectory gives 
one turn after q iterates of the billiard map, so the variable ^ 
should be increased by one if s is increased by / = Length(2). 
Therefore, 


C = C 


f-^ds= 

Jo ds Jo Jq 


ds. 


(14) 


Relation (|9]l is obtained by joining equations (fOl) and (fT4l i. Let 
s = so(^) be the inverse of the solution ^ = ^o('S) of the dif¬ 
ferential equation (0 determined, for the sake of definiteness, 
by the initial condition ^o(O) = 0. By abusing the notation, 
let a:(^) - a:(so(^)) be the curvature in the new angular vari¬ 
able ^ 6 R/Z. Then /c(^) is a 1-periodic real analytic function 
which does not vanish on the reals. Let us assume that there 
exists (5 > 0 such that /c(f) is analytic and does not vanish on 
the open complex strip Is and has singularities and/or zeros on 
the boundary of Is- Note that we are avoiding not only singu¬ 
larities but also zeros of the curvature /('(^). On the one hand, 
the results found by Marvizi and Melrose only hold for smooth 
strictly convex curves, so the zeros of the curvature are a source 
of potential problems. On the other hand, several positive and 
negative fractional powers of the curvature appear in the pre¬ 
vious computations (see also below), and such powers are not 
analytic at the zeros of the curvature. 

Following the numerical evidences in the splitting problems 
of the polynomial standard maps, we thought that r = 27rb, but 
our experiments disprove it. We have obtained that r < 2ji6, the 
equality being an infrequent situation. 

An explanation of such discrepancy is the following one. 
Setcro(^) = (.o(^),r„(^)), ro(^) = p^fi^l2 = CK^I\^)l2q- 
We know that the billiard map (fTTI) can be analytically extended 
to (R//Z) X [0,7r); see Q Proposition 5]. However, we do 
not know whether it can be analytically extended to a complex 
neighborhood in (C//Z) x C of croila) as a —> and q -Hoo 
or not. Hence, the inequality r < 2nb does not look so bad 
in the light of the previous discussion about the Fontich-Simo 
Theorem. In fact, it is commonly accepted that the magnitude 
involved in the exponent of the exponentially small formulas 
for splitting problems is not the minimum distance to the real 
line of the set of singularities of the time parametrization of the 
separatrix but the minimum distance to the real line of the set 
of singularities of the perturbation of the system when evaluated 
on the time parametrization of the separatrix- See ll2^ 38 1 for 


some examples. It seems reasonable to think that one has to 
compute the singularities of the Lagrangian evaluated on the 
solution of (l9]l, which, in its turn, reduces to the study of the 
singularities of •y(so(^))- This is a work in progress. 


4. Model tables 


or, equivalently, as 

C-^ = K^l\s). (13) 

di 

We only need the following observations to determine C. We 
are looking at the (1, q')-periodic trajectories inside Q- We have 


We restrict our study to the perturbed ellipses and perturbed 
circles given implicitly in (|2]l. To be precise, the algebraic curve 
x^+y^ Ib^+ey" - 1 has several real connected components when 
n is odd. Henceforth, we only consider the one that tends to the 
ellipse (or circle) x^ + y^lb^ = 1 as e tends to zero. 
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Let e„ — e„(b) be the maximal positive parameter such that 

Q is analytic and strictly convex for all e e (0, e,,)- (15) 

On the one hand, I„ - (0, +oo) when ii is even. In such cases, 
we will reach the value e = 1 in some numerical computations. 
On the other hand, if ii is odd, the algebraic curve defined by 
+ y^/b^ + ey" = 1 has a singular point on the y-axis when 

e = = e„(b) := 2(n - (16) 

Thus, Q is no longer analytic when e = Our computations 
suggest that e„ = e„ so we restrict our experiments to the range 
0 < e < e„. We note that e^ib) ^ 0.3849/fe^ es{b) ^ 0.1859/fe^ 
and en{b) ^ 0.1232/^?^. We also restrict our experiments to the 
degrees 3 < « < 8. 

The symmetries of our model tables simplify the search of 
some periodic trajectories. If n is even, Q is symmetric with 
respect to both axis of coordinates, so Q is bi-axisymmetric. If 
n is odd, Q is symmetric with respect to the y-axis only, so Q 
is axisymmetric but not bi-axisymmetric. We say that a billiard 
trajectory is axisymmetric when its corresponding polygon is 
symmetric with respect to some axis of coordinates. We only 
compute axisymmetric periodic trajectories, APTs for short. 

First, let us focus on the case odd n. The axisymmetric tra¬ 
jectories inside Q are characterized as the ones with an impact 
point on or with a segment perpendicular to the y-axis. The 
APTs are characterized as the ones satisfying twice the former 
condition. Thus, there are four kinds of APTs inside Q. Be¬ 
sides, only two of these kinds are possible depending on the 
(parity of the) period q. 

The classification for even n is richer because the symmetry 
with respect to the x-axis plays the same role. See Table |2] 

We wanted to study the differences and the Mather’s 
AVTi/^, but instead we will compute the signed differences Dq 
between the lengths of the (l,^)-APTs. Clearly, \Dq\ < 

In some cases, all periodic trajectories are axisymmetric, and so 
A(1.?) = hWxiq = \Dq\. See Proposition|2] 

We will fix the semi-minor axis b and the degree n in our 
numerical experiments. That is, we will study the dependence 
of Dq - Dq(e) on the perturbative parameter e and the period q. 
The quantity Dq{e) is analytic at e = 0 because all (l,q')-APTs 
are so. On the contrary, the period q is a singular parameter of 
this problem because Dq is exponentially small in q. Thus, we 
will deal with: 



• The regular case, where we study the asymptotic behav¬ 
ior of Dq(e) when e —> 0 and q > 3 is fixed; and 

• The singular case, where we study the asymptotic behav¬ 
ior of Dq{e) when q -too and e 6 R is fixed. 


Table 2: Classification of (1, ^)-APTs inside bi-axisymmetric and axisymmetric 
billiard tables Q. In each case, the difference Dq is the length of the (1, ^)-APT 
in red minus the length of the (1, ^)-APT in blue. The gray lines denote the axis 
of symmetry. 


We will see that the classical Melnikov method is suitable 
to study the regular case but it is not so to study the singular 
one. Besides, the Melnikov method gives more information on 
perturbed ellipses than on perturbed circles. The singular case 
is only studied numerically. 


5. Perturbed ellipses 

In this section we restrict ourselves to the case 0 < b < 1. 
We begin with the regular case, so the semi-minor axis b, the 
degree n >3, and the period q>3 are fixed, whereas e —> 0^. 
Since the quantity = A*'’'?^(e) is analytic and vanishes at 
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e = 0 , then 


Ad.?) ^ + O(e^), (17) 

for some coefficient 6 R. This coefficient can be com¬ 
puted by using a standard Melnikov method. In fact, the model 
tables (|7]) have been chosen in such a way that the asymptotic 
behavior of can be determined. The analytical results for 
in the cubic and quartic perturbations are stated below, but 
we need to introduce some notation first. 

Given m € [0,1), the complete elliptic integral of the first 
kind is 

pnH 

K^K(m)^ (l-msin^ey^'^dO. (18) 
Jo 

We also write K' - K'{m) = K(1 - m). 

Proposition 2. Ifb e (0,1) and q > 3, the following properties 
hold. 

1 . Aj*’^' = 0 , for odd n and even q. 

2. There exist some constants c, M-i, M 4 , A 4 > 0, depending 
only on b, such that 


A 


(!.?) 

1 


for n — 3 and odd q, 

. Knqer'^^^^', for n — A and odd q, 
Miiqer‘^‘t, for n — A and even q. 


when q -too. Besides, K 4 — 2 M 4 , and 

nK{b^) _ nK'{\ - b^) 
^ ~ 2K{\-y-) ~ 2K{\-b^) ■ 


(19) 


( 20 ) 


3. If n — 3 and q is odd or if n — A, then there exists 
e„ = e„(b,q) 6 /„ such that all {\,q)-periodic billiard 
trajectories inside O are axisymmetric when e 6 ( 0 , e„). 
In particular, = AWi/^ = \Dq\for all e 6 (0, e„). 


See Appendix A for the proof. The explicit values of M 4 
and M 3 can be found in (lA.lOb . Related computations can be 
found in 0 . 


Remark 1. Similar results hold for any degree n > 5, although 
it is more cumbersome to compute the Melnikov constants M„ 
(and Kn if n is even) and the Melnikov powers mn such that 


M„q"’"e-‘^‘i, 

for odd n and odd q. 


for even n and odd q, 

M„q”'"e-‘^f 

for even n and even q. 


as ^ » H-oo. The Melnikov exponent c does not depend on n. 
From the first order formula (fTTl i. we deduce that 
lim = 1 , 


for any fixed q >3. Next, we wonder whether the roles of e and 
q are interchangeable; that is, if 

lim fA 0 '«VeA ^‘'‘'1 = 1 , ( 21 ) 

£7^+00 L ^ -I 


for any fixed but small enough e > 0 . 

We should compute for big periods q in order 

to answer this question, but instead we compute |D^|/eAj*’®\ 
Both quotients coincide if e is small enough, see Proposition|2] 
We do not compute |D^|/eAj*’'*' for the cubic perturbation and 
even periods, because Aj*’^* = 0 for n - 3 and even q. 

We show the results obtained for the cubic and quartic per¬ 
turbations in Figure [3 These figures are obtained by taking the 
semi-minor axis b - 4/5. Other values for the semi-minor axis 
give rise to similar figures. 

The Melnikov method does not predict the asymptotic be¬ 
havior of in the singular case. That is, limit (l 2 Tb does not 
hold. Indeed, if we fix any e > 0, then the quotient |D^|/eAj*’^' 
drifts away from one as q grows. The drift appears earlier for 
odd periods in the case of the quartic perturbation. As e gets 
smaller, the drift appears at larger periods q. Since the com¬ 
puting time grows quickly when q grows, the computations to 
see that drift when e is very small are unfeasible with our re¬ 
sources. This happens, for instance, when n - A and e - 10“^°. 
See Figure[3] 

Based on these numerical experiments, we guess that there 
exist some critical exponents v„ > 0 such that 


Ad.?) ^ . 


M„eq'""e-‘^f 

for odd n and odd q. 


for even n and odd q. 

M„eq"'"e-‘^f 

for even n and even q. 

H-oo, and V > 

v„. Here, M„, Kn, mn, and 


when e - Oiq^y, q 
c are the Melnikov quantities introduced in Proposition |3 and 
Remark [T] We do not give an asymptotic behavior when n is 
odd and q is even because we do not have any Melnikov predic¬ 
tion for that case. Results about exponentially small asymptotic 
behaviors based on Melnikov predictions are common in the 
literature. For instance, the rapidly forced pendulum is studied 


in 1 ^, 

inJio, 


27 


^ ^ ] and some perturbed McMillan maps are studied 


2811311]. 


Nevertheless, we are interested in a more natural problem. 
Namely, the asymptotic behavior of when q —> - 1-00 and 
e is fixed. As we have said before, we compute the signed 
difference Dq instead of We have numerically checked 

that, if e is small enough, then there exist a constant A 0, a 
power m 6 Z, and an exponent r > 0 such that 


Dq X Aq^fi-'-f 


( 22 ) 


as ^ » H-oo. In fact, the real behavior is slightly more compli¬ 

cated, since these three quantities depend on the parity of q. We 
summarize our results as follows. 


Numerical Result 3. Fix b e (0,1) and n > 3. Let I„ be the 
maximal interval defined in (EH*. There exists e„ — en(b) 6 /„ 
such that the billiard inside 0 verifies the following properties 
for all e € (0, e„). The Borel transform 0 has a radius of con¬ 
vergence p 6 (0, H-oo). Set r — p/2. There exist two constants 
A,B + 0 such that 

{ Bq^^e^^’’'^, for even n and odd q, 

H ^ J ^ 23 ) 

Aq ^e otherwise. 
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Figure 4: The normalized differences Dq tend to a constant when q —» +oo 
in the ranges 1/2 < b < 9/10 and 0 < 6 < 1/10 for the cubic and quartic 
perturbations. If n is even, then we have to study the even and odd periods 
separately. 


Figure 3: The quotient |/)^l/6Aj^’^^ versus the period q for b = 4/5. 
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as q ^ +CX). The quantities p, r, A, and B depend on b, e, 
and n. The constant B is defined only when n is even. Besides, 
limf_>o r — c, where c is the Melnikov exponent defined in (ESll. 

We stated in Conjecture[T]that the function A(^) that appears 
in the exponentially small asymptotic formula ® is constant 
when the billiard table belongs to a certain open set of the space 
of axisymmetric algebraic curves. Thus, the previous numerical 
result fits perfectly into the conjecture. 

It is interesting to compare the Melnikov formulas ( fT9l l with 
the asymptotic formulas (|2^ . The asymptotic behavior of Dq 
does not depend on the parity of q when n is odd. The expo¬ 
nents c and r play the same role. Finally, the factors q^^ and q^^ 
in (|2^ can not be directly guessed from the Melnikov formulas. 

Let us describe our numerical experiments. First, once the 
exponent r is determined (see next paragraph), we compute the 
normalized differences 

^ ( q^ef''‘^Dq, for even n and odd q, 

yq^s''‘>Dq, otherwise. 

We have checked that these normalized differences Dq tend to 
some constant as o' —> -i-oo in the ranges 1/2 < h < 9/10 and 
0<e<l/10. Figure |4] shows that behavior on three different 
scenarios for b - 9/10 and e — 1/10. 

Let us explain how to compute the exponent r - r(b, e, n). 
First, we assume that the exponentially small asymptotic for¬ 
mula (l 22 Ti can be refined as 

j>0 




for some asymptotic coefficients dj 6 R with do - A 0. This 
assumption is based on similar refined asym ptot ic formulas for 
the splitting of separatrices of analytic maps |22, By taking 

logarithms, we find the asymptotic expansion 


i log {q-'”Dq) ^-r+Uog\y^\^-r+Y, 

^ ^ \j>Q ^ J j>Q 


qlj+l ■ 


for some coefficients aj 6 R. Therefore, we can compute r 
by using a Neville extrapolation method from a sequence of 
differences Dq. The longer the sequence, the more correct digits 
in r. We obtain 15 correct digits with the following choices. We 
fix the perturbed ellipse Q, that is, we fix h 6 (0,1), e e R, and 
n > 3. Second, we fix the class of periods q, so that we are 
on one of the cases of Table |2] That is, q = q{k) - 2k + \, 
q - q{k) - 4k + 2, q - q{k) - 4k, ox q - q(k) — 2k. Then, 
we compute Dq with at least 400 correct digits on an increasing 
sequence of 500 periods qi - ^(k,), withk, = ko-i-10/. The initial 
period qo is chosen to be big enough so that \Dqf\ < 10-3000 Jjj 
fact, we perform the Neville extrapolation with two different 
sequences of 500 periods each which allows us to determine 
the number of correct digits in the final result. The power m 6 
{-2, -3) is found by trial-and-error. 

In Figure |5] we display the exponent r - r{e) for several 
values of b for the cubic and quartic perturbations. We also 


Figure 5: The exponent r versus the perturbative parameter e. We also display 
the points (0, c) in solid circles, where c is the Melnikov exponent. We note 
that lime^o+ r = c. Red: b = 1/2. Green: b = 3/5. Blue: b = 7/10. Magenta: 
b = 4/5. Black: b = 9/10. 


depict the Melnikov exponent c at e = 0 in full circles. Note 
that limf^o r - c and r is decreasing in e. 

Next, let us relate the exponent r with the radius of conver¬ 
gence p of the Borel transform (l4]i. Once fixed b 6 (0,1), e 6 R, 
and « > 3, we compute p = p{b, e, n) as follows. 

First, we compute the length of one of the (1, q')-APTs 
inside Q for the same sequences of periods (q,) used for com¬ 
puting Dq. We use a precision of 3000 correct digits in these 
computations. The choice of the APT does not matter, since 
\Dq} < 10“^™ for any period qi > qo. Second, we obtain the 
first asymptotic coefficients Ij in the expansion ([TJ by using the 
Neville extrapolation method again. Third, we determine the 
number of correct digits in each coefficient Ij by comparing the 
results obtained with two different sequences of periods. The 
number of correct digits in Ij decreases as j grows. We always 
get at least 1500 correct digits in /q and at least 40 correct digits 
in / 450 . 

It turns out that the coefficients Ij increase at a factorial rate, 
so the asymptotic series ([TJ is Gevrey -1 and diverges for any q. 
Indeed, we have found that there exist a radius of convergence 


10 




















Figure 6: versus j for b = 9/10, e = 1/20, and n = 4. The dashed 

line corresponds to the limit value p obtained by extrapolation. 


Figure 7: The exponent r (continuous lines with points) and the quantity 27:6 
(continuous lines) versus 6 for b = 4/5. Red: n = 3. Blue: n = 4. 


p - p(b, e,n) > 0 and a constant y = y(b, e,n) > 0 such that 

h " J +“’> 

provided e is small enough. That is, the Borel transform (|4]i has 
a singularity at z = p. In particular, 

p = lim 

j^oo I J " I 

We see this asymptotic behavior in Figure 
The rough approximation 

p ~ I/ 449 /Z 450 I 

only gives about 3 correct digits. If we use an extrapolation 
method based on the asymptotic expansion 

/>o 

the radius of convergence is improved up to 8 correct digits. 
This is the limit value plotted in Figure | 6 ] We stress that this 
asymptotic expansion in powers of is probably wrong since 
the extrapolation becomes unstable after a few steps. 

The radius of convergence p does not depend on the parity 
of the periods of the sequence (qi). Thus, the value of p obtained 
by sequences of different parities must coincide. This provides 
another validation to the number of correct digits of p. 


Proposition 4. Let b e (0,1). Let k(s) be the curvature of the 
unperturbed ellipse E — {(.x,y) e jb^ — \] in some 

arc-length parameter s. Let f 6 be the angular variable 
defined by ( 0 - Let 6 be the distance of the set of singularities 
and zeros of the cur\’ature to the real axis. Then 2n6 — c, 
where c is the Melnikov exponent defined in 


This proposition is proved in Appendix B 


We have numerically checked that the inequality r < 2 tx6 
holds in the ranges 1/2 < < 9/10 and 0 < e < 1/10 for the 

cubic and quartic perturbations. The case b - 4/5 is displayed 
in Figure |7] 


Remark 3. The distance 5 is numerically computed as follows. 
First, we write the curvature k and the length element ds of the 
perturbed ellipse (I7]i in terms of the vertical coordinate y. It 
turns out that there exist three polynomials r(y), p(y) and q(y) 
such that 

^«d, = s(,)dr:=-£3g=d,. 

sjr(y)q(y) 

For instance, r(y) = 1 - y^ jb^ - ey" and deg[p] = deg[g'] = 
2n — 2. Let y± be the roots of r(y) that tend to ±b when e —» 0. 
The points (0,y±) are the vertices on the vertical axis of the 
perturbed ellipse Q. Then 6 - |3f*|/C, where 

C = r ds = 2 r g(y) dy, = f dy, 

Jq Jo 


Remark 2. Taking into account relation r - p/2, we have two 
different ways of computing the exponent r, the direct method 
and the Borel one. The Borel method is computationally much 
cheaper. Indeed, the precision required to compute the differ¬ 
ences Dq. increases along the periods qi whereas it is fixed when 
computing the lengths 

At this point, we have established the relations among the 
Melnikov exponent c, the exponent r, and the radius of conver¬ 
gence p. Next, we relate c with the distance 5 provided by our 
candidate for limit problem, since we are only able to analyti¬ 
cally compute 5 for unperturbed ellipses. 


andy* y± is the root of p(y), q(y), or r(y) that gives the closest 
singularity 6 ’Cj'L to the real axis. That is, y* minimizes 6 . 
The path from y = 0 to y = y* in the second integral should be 
contained in an open simply connected subset of the complex 
plane where the function g(y) is analytic. See Appendix C for 
more details about the function g(y) and their domain of analyt- 
icity, although that appendix deals with perturbed circles only. 


We note that the cusp that appears in the graph of 27rd for 
the quartic perturbation correspond to a perturbative parameter 
e for which two different roots of p(y), q(y), or r(y) give rise to 
the same 6 - |3f*|/C. 
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6 . Perturbed circles 


In this section, we take b — \ 'm the model tables (|2]l. This 
setting is harder than the one of the perturbed ellipses both in 
the regular and singular cases. Let us explain it. 

We begin with the regular case, so we fix the degree « > 3 
and the period q>3 whereas e tends to zero. First, we note that 
the Melnikov exponent c in (l20l i tends to infinity as b tends to 
one, since K{Q) - njl and lim,„^i- K{rn) - +oo. This suggests 
that the Melnikov method gives little information for perturbed 
circles. In fact, in 1331], it is proved that the first order coefficient 
in (fTTT i vanishes for every period q i Q„, where 


Numerical Result 6 . Fix n > 3. Let I„ be the maximal interval 
defined in m- If e 6 /„, then the Borel transform lH has a 
radius of convergence p 6 (0, +oo) . Set r — p/2. There exist 
two non-zero quasiperiodic functions A(q) and B(q) such that 

^ I B(q)q^'^er^''‘i, for even n and odd q, 

‘‘ \A{q'}q~^&^''‘‘, otherwise, 

as q ^ + 00 . Besides, there exists Xn £ K such that 

r-^-^^+Xn+oH) (27) 

n 


f{3,5,...,« - 2,«), for odd «, 

|{2,4,..., n - 2, n) U {2,3,..., n/2), for even n. 


We might use a higher order Melnikov method to look for an 
order k - k{n, q) eN such that 


with 0. This Melnikov computation is not easy so we 

have performed a numerical study instead. As before, we do 
not study but the difference Dg. 

Numerical Result 5. Set 


k — kin, q) — 


2 \ql2n\ 

\2qln\ 

\qln\ 


for odd n and odd q, 
for odd n and even q, 
for even n and odd q, 
for even n and even q. 


Ifn > 3 and q > 2, then there exists dk — djfn, q) 0 such that 
Dgie) = dkf^ + (25) 


This numerical result has two nice consequences on the 
breakup of the resonant caustics of the circular billiard under 
the perturbation + ey" - 1 with any fixed degree n > 3. 

First, all il,q)-resonant caustics break up, because, once fixed 
the period q > 2, 0 for e small enough. Second, there 

are breakups of any order, because the map q kin, ^) 6 N is 
exhaustive. 

We numerically compute the order k in (l25l l by noting that 


k ^ log 


D,ie) \ 

Dqiel&)]' 


as e ^ 0. Finally, there exist a partition /„ = ^ Pn ^ Rn 

satisfying the following properties. 

1. C„ and P„ are open subsets of /„, whereas R„ is a set of 
isolated perturbative parameters. 

2. If e € Cn, both functions Aiq) and Biq) are constant. 

3. If e € Pn, both functions Aiq) and Biq) are periodic. 
Namely, they have the form 

Aiq) - a cosi27T/3q), Biq) -b-\-b co^iAnfq), 

for some average b ^ 0, some amplitudes a,b > 0, and 
some “shared” frequency f > Q. We note that b + fc/2. 

All these numerical results strongly support Conjecture [T] 
For instance, we conjectured that the function Aicf) is either 
constant; Aicf) = a/2, or periodic: Aicf) - acosi2n/3q) in open 
sets of the space of axisymmetric algebraic curves, whereas all 
other cases are phenomena of co-dimension one. This claim 
agrees with the fact that C„ and P„ are open subsets of /,„ 
whereas R„ only contains the perturbative parameters where a 
transition between constant and periodic cases takes place. 

The functions Aiq) and Biq) and the exponent r depend on 
the degree n and the perturbative parameter e, although Biq) is 
defined only for even n. Both functions Aiq) and Biq) “share” 
the frequency in the periodic case. To be precise, the frequency 
of Biq) is twice the frequency of Aiq). It makes sense because 
the exponent in the asymptotic formula containing the function 
Biq) is also twice the exponent in the one containing Aiq). 

The logarithmic behavior of the exponent r stated in (l27l i is 
closely related to the asymptotic formula (l26l) . Indeed, if we 
roughly try to fit the regular behavior (l25l l when e —» 0 with the 
singular behavior = Oiq"'er''‘‘) when q —> -i-oo, then we get 


For instance, \f n - 1, q - 36, and e - 10 then we obtain 
the approximation 

k ^ 5.99999999999999999401..., 


so k - 6. We have tested the formulas listed in Numerical 
Result|5]for all degrees 3 < n < 8 and all periods 3 < q < 100. 
Note that, once fixed n. 


k — kin, q) : 


\ 2 qln, for even « and odd O', 


\q/n, otherwise, 
as O' ^ H-oo. Next, we focus on the singular case. 


(26) 


0(^'"e-"«) = D, = Oi^) - Oie‘>''‘) = 0(e-«''°®"'/"), 


so we guess that r ^ j log e\/n. This reasoning is informal but it 
is confirmed by our experiments. Let us describe them. 

We have set e 6 /„ n Q in all the experiments. First, we 
do so because our multiple-precision computations become a 
bit faster for rational perturbative parameters. There is a sec¬ 
ond reason for that choice. Namely, we change the precision 
very often along our computations, and rational values of e are 
not affected by such changes, because they are stored as exact 
numbers. We have also tried to deal with “big” perturbations in 
order to stress that our results are not perturbative, but we recall 
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(a) Dq versus q for n = 4, 6 = 1, and odd q. 
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(b) Dq versus q for n = 4, 6 = 1, and even q. 



Figure 8: Examples with a constant asymptotic behavior of the normalized differences Dq. 





(a) Dq versus q for « = 3 and e = 1/3. (b) DFT of Dq for « = 3 and 6 = 1/3. (c) Dq - A{q) versus q for « = 3 and e = 1/3. 
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(d) Dq versus q for n = 6, e = 1, and even q. 



(e) DFT of Dq for n = 6, e = 1, and even q. 
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(f) Dq - A{q) versus q for n = 6, e = 1, and even q. 





(h) DFT of Dq for n = 6, 6 = 1, and odd q. (i) Dq - B{q) versus q for n = 6, 6 = 1, and odd q. 


Figure 9: Examples with a periodic asymptotic behavior of the normalized differences Dq. We recall that A{q) = aco?,{2n^q) and B{q) = b + b cos{Anpq). Besides, 
a « 29.4849 and 1/8 in Figure [9^ a ^ 53.2369 and (5 « 0.04614 in Figure [9^ and b ^ -4.9257, b « 7.80853, and jS « 0.04614 in Figure [9^ 
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that e should be smaller than the singular value (fThl l when n is 
odd. 

First, we compute the exponent r - pjl by using the Borel 
method, since it is computationally cheaper than the direct one. 
See Remark |2] Besides, it is not clear how to adapt the direct 
method when the functions A{q) and B{q) oscillate. We follow 
the same steps as in the case of perturbed ellipses. However, the 
Neville extrapolation is more unstable for perturbed circles. In 
order to overcome this instability, now we take sequences {qi) 
of 1000 periods such that 

Once we find r, we compute the normalized differences 
already introduced in (l24li . We have checked that there exist 
two non-zero quasiperiodic functions A{q) and B{q) such that 



for even n and odd q, 
otherwise. 


as ^ » H-oo. 

Some paradigmatic examples of the asymptotic behavior of 
the normalized differences Dq are displayed in Figures |8] and |3 
All these examples are generic in the sense that a small change 
of the perturbative parameter e does not produce any qualitative 
change in the pictures. 

For instance, we see three examples where Dq tends to some 
constant as o' —> H-oo in Figured The constant is A in the second 
and third subfigures, and B in the first one. 

We display a first example of periodic asymptotic behavior 
in Figure [9(a)| for the cubic perturbation and e — 1/3. This 
value e - 1/3 is relatively close to the value 63 ( 1 ) ~ 0.3849 
where the algebraic curve -H -H ey^ - 1 becomes singular. 
Next, we compute the discrete Fourier transform (DFT) of the 
last terms of the sequence Dq. To be precise, the terms in the 
range 10000 < q < 12000 for n - 6 and even q, and in the range 
5000 < q < 6000 otherwise. We discard the first terms because 
Dq “ A{q) and Dq ^ B{q), so the last normalized differences 
are closer to the periodic functions we want to determine. 

The DFT of the normalized differences Dq suggests that the 
periodic function A{q) has a dominant harmonic with amplitude 
a 29.4849 and frequency /? a; 0.375 = 1/8 when e = 1/3 
and « = 3. See Figure [9(b)| This explains why we see eight 
waves in Figure [9(a)l each one with frequency \/3 - 1/8|. This 
situation is a source of problems for the following reason. Let 
us assume that, due to time or computational restrictions, we 
only compute the normalized differences for periods of the form 
qi - qQ + 8i. In that case, we would only see one wave and we 
would get a wrong frequency. The moral of this story is that 
we have to compute the normalized differences for all periods. 
Then we compare the normalized differences Dq with the cosine 
wave A(q) - a cos{27:/3q) as q ^ h-cm. The amplitude a and 
the frequency /? are determined by mixing several tools: the 
DFT, some direct algebraic computations, etcetera. The plot in 
Figure [9(^shows that 


are the following ones. First, we confirm that the frequency of 
the periodic function B{q) is twice the frequency of the cosine 
wave A(q). See Figures |9(e)| and |9(h)| Second, the average of 
B{q) is not zero. This is a surprise, because both the periodic 
functions obtained in similar splitting problems and the periodic 
function A{q) obtained in this billiard problem have generically 
zero average. Third, B{q) - b+bcos{4nPq), but b + bjl, which 
sets another difference with the known asymptotic behaviors for 
splitting problems. 

Next, we present some results about the transition between 
the two generic —“constant” and “periodic”— asymptotic be¬ 
haviors of the normalized differences Dq. That is, we intend to 
visualize what happens at some e* e dC„ n dP„ c R„. 

We focus our attention on the sixtic perturbation: n - 6 . 
Then the normalized differences have “constant” and “periodic” 
asymptotic behaviors for e - 1/10 and e = 1, respectively. We 
study the quantities Dq in a fine grid of perturbative parameters 
in the interval [1/10,1]. Both functions A{q) and B{q) change 
at the same transition value e». Indeed, 


[1/10,23/200] c Ce, [3/25,1] c Pg, 


so the transition takes place at some e* 6 (23/200,3/25). 

Unfortunately, a more precise computation of e* is beyond 
our current abilities, because we do not have a limit problem 
whose complex singularities allow us to determine analytically 
the transition values. An example of such a naly tical computa¬ 


tions for splitting problems can be found in Il23 


Therefore, we only display the normalized differences Dq 
for e - 23/200 and e — 3/25 in Figures [TOl and fTTI to see the 
transition of the functions A{q) and B{q), respectively. 

Let us present some numerical results about the logarithmic 
growth (IZTl i of the exponent r. We have computed the exponent 
r - p/2 by using the Borel method for 3 < n < 8 in a sequence 
of perturbative parameters of the form e, = 2 ^^ jlQ with j > 0 . 
We have plotted the results in Figure [12] On the one hand, the 
curves in Figure [T 2 (a)| look like straight lines with slopes 1 /n, 
as expected. On the other hand, the curves in Figure [T2(b)| tend 
to some constant values Xn > 0. This ends the numerical study 
of such phenomenon. 

Finally, we see that our candidate for limit problem captures 
this logarithmic behavior, although it may not give the exact 
value of the exponent r. 


Proposition 7. Let n > 3 and e € /„. Let k{s) be the curvature 
of the strictly convex curve Q — ((y,y) 6 -H -H ey" — 1) 

in some arc-length parameter s. Let f 6 M.j'L be the angular 
variable defined by (|9|l. Let 6 be the distance of the set of sin¬ 
gularities and zeros of the curvature K{f) to the real axis. There 
exists 6 R such that 


2n5 - 


llog^l 

n 


+ q„+ 0 (e^^" log e). 


(28) 


lim (Dq - A{q)] - 0. 

^^+00 ' ' 

We study the case n - 6 and e = 1 in Figures |9(d)f|9(i)] 
The most interesting phenomena shown up by those pictures 


0 +. 


The proof of this proposition is placed in Appendix C 


The Constantin in (i27l i is always smaller than (or equal to) 
the constant q„ in (i28l l. We compare both constants in Tabled 
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(a) 6 = 23/200. (b) 6 = 3/25. 


Figure 10: Transition of the function A{q) from constant to periodic. We plot the normalized dilferences Dq versus q for /? = 6 and even periods. 



Figure 11: Transition of the function B{q) from constant to periodic. We plot the normalized dilferences Dq versus q for n = 6 and odd periods. 




(b) r — I log 6|/n versus [ log 6|. 

Figure 12: Logarithmic growth of the exponent ?■ as 6 —^ O'*'. Red: n = 3. Green: n = 4. Blue: n = 5. Magenta: n = 6. Cyan: n = 1. Black: n = 8. 
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n 

Xn 

Jin 

3 

0.30... 

1.1358418243... 

4 

0.17... 

1.0703321545... 

5 

0.15... 

0.1488295936... 

6 

0.15... 

1.0385641059... 

7 

0.18... 

0.1823551667... 

8 

0.19... 

1.0332248276... 


Table 3: The constants and r]„, with^„ < rjn, that appear in formulas (H) 
and j28t . respectively. 



Figure 13: The exponent r (dashed lines with points) and the quantity Ind 
(continuous lines) versus |log6|. Red: n = 3. Green: n = 4. Blue: n = 5. 
Magenta: n = 6. Cyan: n = 1. Black: n = 8. 


Constants Xn are computed from the numerical data used 
in Figure |12(b)[ Constants rjn are computed by using the tech¬ 
niques explained in Remark[3] On the one hand, we obtain just 
two significant digits for the constants . On the other hand, 
we can compute rj,, with a much higher precision; here we have 
just written their first ten decimal digits. We see that Xn < T]n 
for n 6 {3,4,6,8). We do not discard the equalities ^5 = 775 and 
Xi = rjj. In order to elucidate them, we compare the exponent 
r with the quantity Ind, as we have done before for perturbed 
ellipses at the end of Section |5] The results are displayed in 
Figure [T3] where we see that our candidate for limit problem 
gives the exact exponent r in two cases. 

To be precise, our numerical results suggest that: 

• If n 6 {3,4,6,8), then r < 2 tx6 for all e 6 (0,1/10); and 

• If n e {5,7}, then r = Ind for all e 6 (0,1 /10). 
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Appendix A. Proof of Proposition |2] 


We will use many properties of elliptic functions listed in 
the books 1142114311 . a couple of technical results about ellip¬ 
tic billiards contained in i^l, and the subharmonic Mel¬ 
nikov potential of billiards inside perturbed ellipses introduced 

inllsl. 

We consider the unperturbed ellipse 

£ = {(x,y) 6 = 1 }, 0<h<a. (A.l) 

It is known that the convex caustics of the billiard inside E are 
the confocal ellipses 


Ca = ux,y) 6 


p2 . 




2 ^ — — A^ 


= 1 


0 < A < b. 


There is a unique (p, q')-resonant elliptic caustic Ca for any rela¬ 
tively prime integers p and q such that 1 < p < qjl. The caustic 
parameter of the (p, q')-resonant elliptic caustic is implicitly de¬ 
termined by means of equation (IA.3b . 

The complete elliptic integral of the first kind is 


K = K(m) = 


1 ' 


(1—msin^f) 


Its argument m e (0,1) is called the parameter. We also write 
K' - K'(m) - K{1 - m). The amplitude function ip - amf is 
defined through the inversion of the integral 


t - 


pp 

( 1 - 

Jo 


m sin^ (p) d(p. 


The elliptic sine and elliptic cosine associated to the parameter 
m 6 (0,1) are defined by the trigonometric relations 

sn f = sn(f, m) = sin tp, cn f = cn(f, m) = cos tp. 

If the angular variable (p changes by 2n, the angular variable t 
changes by AK. Thus, any 27r-periodic function in ip, becomes 
4£-periodic in t. By abuse of notation, we will also denote 
the 4A-periodic functions with the name of the corresponding 
27r-periodic ones. For example, if q{ip) - {a cos ip, b strap) is 
the natural 27r-periodic parameterization of the ellipse E, then 
q{t) - (a cn t, b sn t) is the corresponding 4A'-periodic parame¬ 
terization. The billiard dynamics associated to an elliptic caus¬ 
tic Ca becomes a rigid rotation f i-» f H- b in the variable t. It 
suffices to find the shift 6 and the parameter m associated to each 
elliptic caustic Ca- The parameter m is g iven in 0441 Eq. (3.28)] 
and the constant shift b is given in 0441 p. 1543]. We list the 
formulas in the following lemma. 

Lemma 8. Once fixed an elliptic caustic Ca with A 6 (0, b), the 
parameter m 6 (0,1) and the shift 6 6 (0, 2K) are 

„2 u 2 p »/2 

b/2= {l-msin^f)-^/^d(l>, (A.2) 


2 ^ — A^’ 


where & 6 (0,7r) is the angle such that sm(d/2) — Ajb. The 
segment joining the points q(t) and q(t + b) is tangent to Ca for 
all f 6 R. 
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From now on, m and S will denote the parameter and the 
constant shift defined in (IA.2b . Observe that the elliptic caustic 
Ca is ip, q')-resonant if and only if 

qd = 4-Kp. (A.3) 

This identity has the following geometric interpretation. When 
a billiard trajectory makes one turn around Ca, the old angular 
variable ip changes by 2 n, so the new angular variable t changes 
by 4K. Besides, we have seen that the variable t changes by 
b when a billiard trajectory bounces once. Hence, a billiard 
trajectory inscribed in E and circumscribed around Ca makes 
exactly p turns after q bounces if and only if (IA.3b holds. 

We consider the elliptic coordinates {p, ip) associated to the 
semi-lengths 0 < b < a. That is, (p, p) are defined by relations 

X - cr cosh p cos p, y - cr sinh p sin (A.4) 

where cr = is the semi-focal distance of E. The ellipse 

E in these coordinates reads as p = po, where cosh/ro = ajcr 
and sinhjUo = bjcr. Hence, any smooth perturbation of E can 
be written in elliptic coordinates as 

p^ Pa+ epi ip) + 0(e“), (A.5) 

for some 2;7r-periodic function p\ : R ^ R. 

Lemma 9. Let p and q be two relatively prime integers such 
that I < p < ql2. Let Ca be the {p,q)-resonant elliptic caustic 
of the ellipse ( L4.7b . Let 


Proof. It follows directly from (jr^, Proposition 10]. □ 

Lemma 11. The following properties hold for p\{ip) — cos^ p. 

1. The Melnikov potential L^^’‘‘\t) has just two real critical 
points (modulo its periodicity), none of them degenerate. 

2. There exist an exponent Ij — a, b) > 0 and a quan¬ 
tity Q4 = Q4(a»t, a, b,p,q) > 0 such that 

* lQ4e“^^, for even q, 


as piq (JL>„ e {0) U ((0,1) \ Q). 

3. There exist r4 = Yi\(a)t,a,b) > 0 and ©4 = & 4 (a,b) > 0 
such that 


Q. 4 {tOt,a,b,p,q) 


r4q^, if CJ, 6 (0,1) \ Q, 
&4Pq, ift^t - 0. 


4. ^(0,a,b) = TtK'il - (bla)^)l2K(l - (b/af). 
Proof. By definition, if p\{p) - cos^ p, then 


Lf^V) = 2T^cn2(r + ;b). 


>=o 


The square of the elliptic cosine is an elliptic function of order 
two, periods 2K and 2K''\, and double poles in the set 


P^K'i-v 2X1. + 2K'iL. 


be the maximal difference among lengths of(p, q)-periodic tra¬ 
jectories inside the perturbed ellipse (O- Let Pi (t) be the 4-K- 
periodic function associated to the 27t-periodic one p\{p). Let 

L^l^'‘^\t)^2AY^pi(t + j 6 ). (A.6) 

J=o 


be the subharmonic Melnikov potential of the caustic Ca for 
the perturbed ellipse (E3D- If L^‘^\t) does not have degenerate 
critical points and e > 0 is small enough, then there is a one-to- 
one correspondence between the critical points of L^f’‘^\t) and 
the (p, q)-periodic billiard trajectories inside (O- Besides, 


^ max-min 


Proof. It follows directly from results contained in Il34ll . □ 


We will determine the asymptotic behavior of First, 

we study the asymptotic behavior of the (p, ^)-resonant caustic 
Ca as pjq —> 0^. 


Lemma 10. If Ca is the (p,q)-resonant elliptic caustic of the 
ellipse (EZl), then A ^ ’L.pjq as pjq —» 0^, where 


E - E(a,b) ab J" ^s(s - b^)(a^ - s)^ ^ ds. (A.7) 


Besides, the principal part of any pole t e P is —nr^{t — t)~^. 
In particular, l!f’‘^\t) is also an elliptic function of order two, 
and so, it can be determined (modulo an additive constant) by 
its periods, poles, and principal parts. 

We study the cases odd q and even q separately. 

If q is odd, then l!'f’‘^\t) has periods 2Klq and 2K'i and 
double poles with principal parts -2/lm '(f - in the set 

'IK 

P = K'i + —Z -H 2K'i. 

q 

It is known that K'{m)IK(m) is a decreasing function such that 


,. K'(m) 

lim -= -Hoo, 

m-»o+ K(m) 


K'im) 

hm - = 0. 

m->i- K{m) 


Therefore, there exists a unique m^ e (0,1) such that 


K'g K'iniq) K'(m) K' 

Kq K{mq) ~ ^ Kim) K' 

Henceforth, we write that K - Kim), K' - K'im), Kq - Kiniq), 
and K'q = K’iniq) for short. Thus, 

= const. -V 2AiqKqlK)^imqlm) cn^iqKqt jK,mq), 


which has just two real critical points (modulo its periodicity), 
none of them degenerate. Besides 

= maxL^’''^ - min if = 2 AiqKqlK)\mqlm). 
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If piq ait e (0,1) \ Q, then q +cx) and d ^ d* 6 {0,b), 
where Ca, is the elliptic caustic with rotation number w*, so 


a^-b^ 

nit ^ 6 (0,1), 


K' 

~K 


Kt 


a? — At 

K'init) 


K{mt) 


6 (0, +oo), 




0+, 

7T 

2 ' 


Using 17.3.14 & 17.3.16], we get the asymptotic formula 


16e where ^ nK^jlKt. Finally, we obtain that 






nitK: 


^q^e as pjq 


Ut and q is odd. (A. 8) 


If q is even, then cn^(f + q6l2, m) - cn^(f, m) and 

q/2-l 

Lf’‘>\t) = 4d ^ cn^(f + j6, m). 


so L 


j=o 

has periods 4K/q and 2fif'i. In this case, 
4n^At 


i(P,9) 

^1 


‘-q^e as pjq ai* and q is ( 


(A.9) 


Next, we study the case w* = 0, when the {p, q')-periodic 
orbits approach the boundary. In this case, 

, A 1 / n2 >- >-.a nK'{\-{blaf) 

At^O, nit^l-ib/a), f = ^(0, a,/?) = —-— . 

2 a (1 - (blay) 

Since At - 0, we need the asymptotic behavior of the caustic 
parameter d as pIq ^ 0^. We recall that d x 'B.plq in that case, 
where S = H(a, b) is the integral defined in (IA.7b . Hence, 

47r^dt ^ 47r^H(fl, b) 


r. = 


04 = 


nitKi ’ (1 - (b/a)^)K( 1 - {b/a^ ’ 

and this ends the proof of the lemma. □ 

Lemma 12. The following properties hold for pi(ip) — — sim^. 

1. Ifq is even, then = 0 ‘^nd = 0. 

2. If q is odd, then has just two real critical points 

(modulo its periodicity), none of them degenerate. 

3. Let f(aJt, a, b) be the exponent introduced in Lemma m 
Ifq is even, then there exists — Q.i(a>t,a,b, p,q) > 0 
such that 


V (p.q) 


113 e 




pIq -> ai» e {0) U ((0,1) \ ' 


4. There exist r3 = T 2 (a>t,a,b) > 0 and 03 = ®^{a,b) > 0 
such that 


Q.i(tOt,b,a, p, q) - 


\T^q, i/w. 6 (0,1) \ I 

l03P, ifa>t-0. 


Proof If q is even, then p is odd, sn(f + b/2) = -snf, and 
Lf^^(f) = -2di:p;sn(f + ;b) = 0. 

The case odd q follows the lines of the proof of LemmafTTI 
The constants are 

87rd» „ 87rS(a, b) 


r3 = 


sJinlKt ’ 


03 = 


(l-(/7/a)2)i/2/:(l_(/,/a)2)’ 


where Ca, is the elliptic caustic with rotation number w*, m* = 
(a^ — b^)l(a^ - d,), and Kt - K(mt). We omit the details. □ 


Next, we relate the original perturbed ellipses (|7]i written in 
Cartesian coordinates, to the perturbed ellipses (lA.Sl l written in 
elliptic coordinates. 

Lemma 13. Set 0 < b < a. 

1. The perturbed ellipse ( L4.5I) with pi{ip) — — sin^ has, up 
to terms of second order in e, the implicit equation 


x2 (y-eb^/a)^ - b^ 3 , 

b^ 


2. The perturbed ellipse (O with p\{ip) — cosf tp has, up 
to terms of second order in e, the implicit equation 

^2 2 2_^2 

for some semi-lengths a — a + 0(e) and fi - b + 0(e). 

Proof It is a tedious, but straightforward, computation. □ 

Finally, we get the claims stated in Proposition |2] from the 
previous results by using that a - a + 0(e) and (3 - b + 0(e) 
and by taking a = 1. To be precise, then 


c = c(b) = ^(0, l,b) = 

M3 = Miib) = 

M4 = M^(b) = 


nK'{\-b-) 

2K{\ - b^) ’ 
b''®A,(\,b) 47Tb^'R{\,b) 


2{\-b^) (1 -/72)3/2^(1 _^2) 

b^®A{l,b) 2n^b^'E.(\,b) 


, (A. 10) 


2(1-172) (1 -l72)2/:(l -/72)2’ 

where the elliptic integral S = S(a, b) is defined in (IA.7b . 

Appendix B. Proof of Proposition |4] 

We parameterize the ellipse by using the angular variable (f. 
That is, we use the parametrization (j{ip) - (cos (f, b sin ip). The 
curvature of the ellipse E at the point cr{(p) is 


K{ip) = 


1 


(sin^ ip + b'^ COs2 ^)2/2 ^2( J _|. y 5 jjj2 i^p/2 


where v = (1 - b^)lb^ > 0. The arc-length parameter s and the 
angular parameter ip are related by 


di 

dip 


(ip) - ||cr'((^)|| = yjsin^ ip -i- b^ cos2 ip - b -\Jl + v sin^ ip. 


First, we compute the constant 


rtnjl 


C= r ds ^ 4b-^'^ f (\+vti\n^p)-''^dp 
Je Jo 

= 4b-^'^K(-v) = 4b^^^K(l - b^). 


We have used 0421 17.4.17] in the last equality. 

The incomplete elliptic integral of the first kind with ampli¬ 
tude p 6 (0,7r/2) and parameter m 6 (0,1) is 

pp 

F(p\m) - I (1 - OT sin^ 6)“*^^ d0. 

Jo 
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This definition can be extended to complex amplitudes and any 
real parameter ll^ . Note that F{nl2\m) — K{m). 

The curvature K((p) has no complex zeros but has complex 
singularities at the points such that sin^ ^ = -1 /v. This equa¬ 
tion becomes sinh^ if/ - I jv under the change ip — iip. Let be 
the only positive solution of the previous equation. Any singu¬ 
larity of K{ip) has the form 

ip — ipfj ±ii//t + tin, « e Z. 

Let ^ be the complex singularity of k{^) associated to ipf; 
through the change of variables 




Jo Jo 


The complex path in this integral is the segment from 0 to ip. 
Next, we compute the complex singularities : 

e„^c-^ r i^i\0)^{0)d0 

Jo 

= C^'r'^V(ii/r.-H «7r| - v) 

= 2nC-^b-^^^K(-v) + iC-^b-^>^F(n/2\b^) 

= 2nC-^b^l^K{\ - b^) + \C-^b^l^K(b^) 

^nl2 + \C-^b^l^K'{\-b^). 


By symmetry, = ~^^n- We have used formula 0421 17.4.3] 
to compute Fiiif/^ + nn\ — v) , formula 11421 17.4.8] to compute 
F{\^t\ - v), and formula ||43 17.4.15] to compute F(nl2\b~). 

Therefore, the distance S of the set of singularities and zeros 
of the curvature k(^) to the real axis is 

Appendix C. Proof of Proposition |7] 

Fix the integer n> 3. We consider the perturbed circles 

e = {(x,y) 6R2 : + = l} (C.l) 

where 0 < e <sc 1 is a small perturbative parameter. 

Let C = C(e) be the constant defined in (|3). If e = 0, then 
2 is a circle of radius one with curvature k = 1, so 

C(0) = I ds = j ds = Length(2) = 27r. 

Jq Jq 

We note that (1C. lb is a smooth perturbation of a circle of radius 
one, so C(e) is smooth at e = 0 and 

C = C(e) = C(0) -H 0(e) = 27r -t 0(e). (C.2) 

We introduce the polynomial r{y) = 1 - - ey". Note 

that (x,y) 6 2 if ^nd only if x^ = r(y). By taking derivatives 


twice with respect to y the implicit relation x^ = r(y), we get 
the auxiliary polynomials 


p(y) 


d^x Ir'iyfY r(y)r"(y) 


dy 

= 1 + epn-iy' 


2 \ 2 
2 

2 


H- epny" + e^p2n-2y^’' 


dx 


q{y)^x + X— =r(y)H- 


dy 

= 1 H- eqny'' + e^qin^ 


/(y) 


,.2n-2 


2y 


whose coefficients are pn-i - «(n-l)/2,p„ = -(n-1)(«-2)/2, 
P 2 n -2 - -n{n - 2)/4, q,, - n- I, and ^2«-2 = n^/4. The length 
element and the curvature at the point (v, j) € 2 are 


ds = 




d^x I, / dx 


dy2 


dy 


2\ 


-3/2 


p(y) 

q^Hy)' 


The curvature should be positive, which explains the minus sign 
in the formula for /c(y). Thus, we can relate any singularity (or 
any zero) y* 6 C of the curvature /c(y), with the corresponding 
singularities (or zeros) 6 C//Z and 6 C/Z by means of 
the formula 


= r ds= f g(y) dy, 

Jo Jo 


where 


g(y) ;= K^'^iyf^iy) 




,2/3(3;) 


Yr{y)q(y) 

Let 7? c C be the union of the complex rays jayo : a > 0], 
where yo is a root of p(y), qiy) or r(y). The function g(y) is 
analytic in C \ 7?, so we will avoid the set 7? when computing 
the integral g(y) dy along complex paths. 

Lemma 14. Let 0 < e 1 and n e N with n>3. 

1. The polynomial p(y) has n roots of the form 

-H 0(e*^"), z" = 2/((n - l){n - 2)); 
and n — 2 roots of the form 
ze”'/(""2) -H 0(e'4""2))^ ^n-2 _ _ 2)ln. 

2. The polynomial qiy) has n roots of the form 

+ Oie^h, z" = -!/(«-!); 
and n — 2 roots of the form 

Zf-i/(«-2) ^ o(e'^<"-2>), z"^2 = -4(n - l)/n\ 

3. The polynomial r(y) has n — 2 roots of the form 

^e->/(«-2) + o(^l/(«-2))^ z"^2^_l. 
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and two real roots of the form y± = +1 -I- 0(e). 















Besides, each one of these roots depends on some positive frac¬ 
tional power of e in an analytic way. 

Proof If wo(z) is a polynomial with a simple root zo and w\{z) 
is another polynomial, then w{z) - wo{z)+pw\{z) has some root 
of the form z - Zo + Oip) which depends analytically on p. The 
roots y± = +1 + 0(e) of the polynomial r(y) - 1 - y^ - ey" are 
obtained directly with wo(z) - 1 - wi(z) = -z", and p — e. 

If we take p - then 

p{e~^'"z) = 1 + p„z" + p{pn- 2 f^^ + Pln-lZ^"^^), 
q{e~^'"z) = 1 + qnz” + pqm-iz^"^^, 

and we find the n roots with an 0(e*^")-modulus of p{y) and the 
n roots with an 0(e^^'')-inodulus of qiy). 

If we take p - e2/(n-2)^ 

pp{e^^l^’'^'^h) = Z"(p„ + pin-iz’'^^) + p{l + Pn- 2 ^^^), 

= z"(^„ + q 2 n- 2 f^^)+P, 

pr(e~^^^"^^h) - -Z^(l +z"^^) +/t, 

and we find the n-2 roots with an 0(e*^^"'^^)-modulus of p{y), 
the n-2 roots with an 0(e*^^"”^^)-modulus of q(y), and the n-2 
roots with an 0(e*^’^"^^*)-modulus of r(y). □ 

Lemma 15. IfO<e<s: 1, neN with n > 3, andyi, e C is a 
root of ply) or qiy) with an 0{e^^^'')-modulus, then there exists 
a constant 77* e R such that 

13^*1 = +77*+0(e2/« log e), 

n 

as e ^ 0^. 

Proof For simplicity, we assume that y* is a root of qiy) such 
that !Ry* < 0 and 3 y* > 0 . Other cases require minor changes. 

If ro = (77 - 1)-'/72, 7-* = e‘/"|y*|, and 0* = argy*, then 
7r/2 < 6>* < nnjin + 1) and r* = 2ro + Oie^^"), because y* = 
e“'^"z + 0(e'^") for some z € C such thatz" = -l/(n - 1) < 0. 

We compute giy) dy by integrating over the path 

cr* = (Ti U 0-2 U 0-3, where 

cTi = {e"'^"if: 0 < f < ro), 

0-2 - {e"'^"roe®‘ ; 7 r /2 < 0 < 0 *), 

0-3 = {e“'^"e®*V : ro < r < r*). 

This path only intersects the set of rays 'R at its endpoint y*, 
since the 2 n roots of piy) and qiy) with an 0(e“'^")-modulus 
have pairwise different arguments when e ^ 0^. 

We write = £* giy) dy = giy) dy = +^2 , where 


- 

1 giy) dy = 

pro 


JtTl 

Jo 


r 

pSi, 

ft = 

1 giy) dy = 

e-'/7oe®‘ig(e-'/"7-oe''‘)d0. 


kj cr2 ' 

J71/2 

ft = 

1 giy) dy = 

r%-i/«e«*i^(g-i/«efl*V)dr. 


Jcr^ 

Jro 


In order to study fi, we consider the function 

hit) e-'/" Vf2 + e2/«^(e-i/nif) ^ /,Q(f) + o(e"^"), (C.3) 


where hoit) = (1 + p„i'’n^^\l + qJ'’t'T^^f The function /70(f) 
is smooth in the interval [0, ro] and hoit) = 1. Besides, 

fi=i f \t^ + e^l’')-^l^hit) df = |i + |i + |i + li, 

Jo 


where 


^1 =i 


‘f 


dt 


= iargsinh(e *^"ro) 


0 Vf^ + e ^ h ' 

iU2i^+ilog(2ro) + 0(e2/''), 


=i 

=i 


‘f 

‘f 


hoit) - 1 


hoit) - 1 


df. 


Vf2 + e2/n 


1 df = 0(e2/”loge). 


■'-o hit) - hoit) 

Vf2 + e2/« 


df = loge). 


The integral is immediate. The integral does not depend 
on e. The integral fi is bounded using ideas from the proof 
of Lemma 23 in OTSll . The integral is bounded using (IC.3l l. 
Hence, we have already seen that there exists 771 e M such that 

|3^i| = ^i^+77i+0(e"/''loge). 


The study of ^2 and ^3 is easier, because 

6 = I2 + 0(e2/”), ft = I3 + 0(e"^"), 


for some constants ft and ft that do not depend on e. 

For instance, ft depends on in an analytic way, because 
the integrand e“*^'Voe®‘ig(e“'^"roe®‘) and the argument 0* are 
both analytic in e^^", and all the singularities of the integrand 
are far from the integration path. The study of ft is similar. □ 
Finally, if b is the distance of the set of singularities and 
zeros of the curvature 7f(ft to the real axis, then 


27rb = — min {|3ft| : y* is a root with an 0(e*^")-modulus| 
^ L!^+^ + 0(e2/"loge), 


where the constant 77 = 77^ e R is equal to the smallest constant 
77* provided by Lemma [15] among all the roots of piy) and qiy) 
with an 0(e*^")-modulus. We have also used relation (1C.2b in 
the last equality. 

We do not care about the roots y± = +1 + 0(e) of r(y), 
since they correspond to points where y is not a true coordinate 
over the perturbed circle Q. To be precise, the points (0, y±) are 
the two vertices of Q over the symmetry line {x - 0), and the 
curvature has a finite positive value at them. Nor do we care 
about the roots whose modulus is because 
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where 


Z’'l^\pn + P2n-2Z’'-^f>^ 

(1 + z"~^y^^{qn + qin^iz'^^^yi^' 

This implies that, if y* 6 C is one of those farther roots, then 

13^*1 = 0 ( 1 )) 

for some constant v* 6 R. That is, the farther roots give rise to 

much bigger imaginary parts. 
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